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Abstract 

This  Special  Report  documents  available  research  pertaining  to  the  disper¬ 
sion  of  bed-load  sediment  transport,  specifically  the  dispersion  coefficient. 
This  knowledge  can  inform  the  evaluation  and  improvement  of  existing 
U.S.  Army  Corps  of  Engineers  (USACE)  sediment  transport  models.  A  dis¬ 
cussion  is  presented  on  normal  and  anomalous  dispersion.  Several  aspects 
of  bed-load  dispersion  are  discussed,  including  whether  the  assumption  of 
normal  dispersion  and  isotropy  are  valid.  Numerical  dispersion  coeffi¬ 
cients  are  relative  to  the  context  of  all  considerations  within  a  modeling 
framework  and  affect  sediment  erosion,  deposition,  and  the  resulting  wa¬ 
ter  quality.  As  witnessed  by  the  number  of  recent  publications,  sediment 
bed-load  dispersion  is  an  active  area  of  research,  and  the  published  values 
for  the  bed-load  dispersion  coefficient  vary  greatly. 

This  document  is  intended  for  scientists  and  engineers  with  a  basic  under¬ 
standing  of  physics  and  numerical  modeling  of  sediment  transport. 


DISCLAIMER:  The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication,  or  promotional  purposes.  Ci¬ 
tation  of  trade  names  does  not  constitute  an  official  endorsement  or  approval  of  the  use  of  such  commercial  products. 
All  product  names  and  trademarks  cited  are  the  property  of  their  respective  owners.  The  findings  of  this  report  are  not  to 
be  construed  as  an  official  Department  of  the  Army  position  unless  so  designated  by  other  authorized  documents. 

DESTROY  THIS  REPORT  WHEN  NO  LONGER  NEEDED.  DO  NOT  RETURN  IT  TO  THE  ORIGINATOR. 


ERDC/CHL  SR-16-5  iii 


Contents 

Abstract . ii 

Figures  and  Tables . iv 

Preface . v 

Unit  Conversion  Factors . vi 

1  Introduction . 1 

1.1  Approach . 2 

1.2  Objective . 3 

2  Background . 4 

2.1  Diffusion  and  dispersion  concepts . 4 

2.2  Dispersion  coefficient . 6 

3  Bed-Load  Dispersion . 9 

3.1  Longitudinal  dispersion . 10 

3.1.1  Normal  dispersion . 10 

3.1.2  Evidence  of  anomalous  dispersion . 14 

3.1.3  Fractional  Advection- Dispersion  Equation  (ADE) . 20 

3.2  Transverse  dispersion . 21 

4  Conclusion . 22 

References . 23 


Appendix:  Definitions  of  Symbols 


30 


ERDC/CHL  SR-16-5  iv 


Figures  and  Tables 

Figures 

Figure  1.  Example  of  2-Dimensional  advection-dispersion  (advection  from  left  to 
right) . 6 

Figure  2.  Example  of  dispersion  that  follows  a  Gaussian  distribution,  while  also 
undergoing  advection  from  left  to  right  (progressing  in  time  from  ti~>  t2~>  t3 ) . 15 

Figure  3.  Example  of  dispersion  that  produces  a  heavy-tail  distribution 

(progressing  in  time  from  ti  ->  t2  t3) . 16 

Tables 

Table  1.  Solute  transport  dispersion  coefficient  equations  found  in  the  literature . 7 

Table  2.  Bed-load  dispersion  coefficient  formulas  and  relations  available  in 

literature . 13 

Table  3.  Bed-load  dispersion  coefficient  values  available  in  literature . 13 


ERDC/CHL  SR-16-5 


v 


Preface 

This  study  was  conducted  by  the  Coastal  Inlets  Research  Program  (CIRP), 
which  is  funded  by  the  Operation  and  Maintenance  (O&M)  Navigation 
business  line  of  the  Headquarters,  U.S.  Army  Corps  of  Engineers 
(HQUSACE)  under  Project#405370,  “Bed-Load  Sediment  Transport  Re¬ 
search.”  The  CIRP  is  administered  for  Headquarters  by  the  U.S.  Army  En¬ 
gineer  Research  and  Development  Center  (ERDC),  Coastal  and  Hydraulics 
Laboratory  (CHL),  Vicksburg,  MS,  under  the  Navigation  Research,  Devel¬ 
opment,  and  Technology  (RD&T)  Portfolio.  At  the  time  of  this  study,  Jef¬ 
frey  A.  McKee,  Chief,  Navigation  Branch,  HQ  was  the  HQUSACE 
Navigation  Business  Line  Manager,  W.  Jeff  Lillycrop,  CHL,  was  the  ERDC 
Technical  Director  for  Civil  Works  and  Navigation  RD&T,  Charles  E.  Wig¬ 
gins  was  the  ERDC  Associate  Technical  Director  for  Navigation,  and  Dr. 
Julie  Rosati,  CHL,  was  the  CIRP  Program  Manager. 

The  valuable  feedback,  reviews,  and  advice  from  Kervi  Ramos,  David 
Abraham,  and  Gary  Brown  are  gratefully  acknowledged. 

The  work  discussed  herein  was  performed  by  the  River  Engineering 
Branch  (CEERD-HF-R)  of  the  Flood  and  Storm  Protection  Division 
(CEERD-HF),  ERDC-CHL;  and  the  Institute  for  Water  Resources,  Hydro- 
logic  Engineering  Center  (IWR-HEC),  Davis,  CA.  At  the  time  of  publica¬ 
tion,  Keith  Flowers  was  Chief,  (CEERD-HFR);  and  Ty  Wamsley  was 
Division  Chief,  (CEERD-HF).  The  Deputy  Director  of  ERDC-CHL  was  Dr. 
Kevin  Barry  and  the  Director  was  Jose  E.  Sanchez. 

The  Commander  of  ERDC  was  COL  Bryan  S.  Green,  and  the  Director  was 
Dr.  Jeffery  P.  Holland. 


ERDC/CHL  SR-16-5 


vi 


Unit  Conversion  Factors 


Multiply 

By 

To  Obtain 

feet 

0.3048 

meters 

inches 

0.0254 

meters 

square  feet 

0.09290304 

square  meters 

square  inches 

6.4516  E-04 

square  meters 

ERDC/CHL  SR-16-5 


1 


1  Introduction 

Research  concerning  coastal  and  riverine  sediment  transport,  specifically 
erosion  and  deposition,  affects  public  safety,  water  resource  management, 
and  operations  and  maintenance  practices.  USACE  developed  and  success¬ 
fully  applied  several  numerical  morphologic  and  sediment  particle  track¬ 
ing  models  (e.g.,  Tate  et  al.  2009;  Sanchez  and  Wu  2011;  Gibson  et  al. 
2010);  however,  existing  representation  of  sediment  transport  is  limited 
and  numerical  calculations  often  vary  by  more  than  an  order  of  magni¬ 
tude.  As  part  of  the  continuing  model  development,  the  present  work  fo¬ 
cuses  on  improving  the  dispersion  coefficient  and  dispersion  equations 
used  for  bed-load  sediment  transport.  For  the  discussion  herein,  sediment 
advection  can  be  understood  as  the  transport  of  sediment  along  with  the 
general  flow  of  the  water,  which  is  typically  related  to  the  average  velocity. 
Sediment  dispersion  describes  how  the  sediment  particles  spread  relative 
to  one  another  as  they  are  being  transported.  Bed-load  dispersion  con¬ 
cerns  the  spreading  of  sediment  particles  that  are  moving  along  the  bed  of 
the  water  body. 

Accurate  simulation  of  the  bed-load  dispersion  is  important  for  under¬ 
standing  and  predicting  the  spread  of  specific  materials.  Applications 
where  users  are  interested  in  tracking  the  transport,  mixing,  and  sedimen¬ 
tation  of  specific  materials,  perhaps  a  specific  grain  size  of  sediment  or  a 
contaminant,  would  benefit  from  improved  dispersion  modeling  capabili¬ 
ties.  For  these  applications,  greater  accuracies  in  estimating  the  arrival 
and  distribution  of  material  concentrations  are  important. 

Historic  sediment  transport  studies  have  focused  on  estimating  the  equi¬ 
librium  sediment  transport  rate  (e.g.,  Yalin  1977;  Meyer-Peter  and  Miiller 
1948)  based  on  solving  the  continuity,  or  Exner,  equation  (Parker  et  al. 
2000)  over  a  control  volume.  Methods,  such  as  Meyer-Peter  and  Miiller, 
computed  mass  change,  but  could  not  capture  the  movement  of  particles 
or  pulses  explicitly.  More  recently,  several  studies  have  focused  on  sedi¬ 
ment  wave  scour  rates  (Abraham  et  al.  2011)  and  particle  velocity  (Lajeu- 
nesse  et  al.  2010;  Furbish  et  al.  2012b),  which  are  the  primary  components 
of  sediment  advective  transport.  However,  there  have  been  relatively  few 
studies  of  sediment  dispersion  (e.g.,  Einstein  1972;  Hubbell  and  Sayer, 
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1964;  Drake  et  al.  1988)  before  the  last  decade,  especially  bed-load  disper¬ 
sion.  Furthermore,  a  general  lack  of  consensus  on  sediment  dispersion 
within  the  literature  exhibits  a  need  for  further  research. 

1.1  Approach 

Recently,  an  ERDC-CHL/IWR-HEC  team  of  Gibson  et  al.  (2016)  con¬ 
ducted  laboratory  experiments  of  bed-load  sediment  transport  to  measure 
dispersion.  They  tracked  colored  sediment  pulses  along  a  flume  under 
steady  equilibrium  sediment  transport.  Bed  samples  were  analyzed  to 
study  bed-load  advection  and  dispersion.  The  experiments  serve  as  a  valu¬ 
able  validation  dataset  for  sediment  transport  models  and  provide  insight 
into  the  physical  processes  of  bed-load  transport.  Two  USACE  sediment 
transport  models,  namely  the  Adaptive  Hydraulics  (AdH)  model  with  the 
sediment  library  (SEDLIB)  and  the  Hydrologic  Engineering  Center’s  River 
Analysis  System  (HEC-RAS),  will  be  tested  with  results  from  Gibson  et  al. 
(2016)  and  the  relationships  discussed  herein,  which  will  inform  any  rec¬ 
ommended  changes  to  the  models. 

AdH  hydrodynamic  model  with  SEDLIB  (Brown  et  al.  2014)  was  devel¬ 
oped  at  the  Engineer  Research  and  Development  Center’s  (ERDC)  Coastal 
and  Hydraulics  Laboratory  (CHL).  In  addition,  the  Adh/SEDLIB  is  under¬ 
going  further  development,  and  it  is  available  to  the  public.  AdH  is  capable 
of  simulating  2-D  or  3-D  hydraulics  along  with  sediment  modeling  and  the 
transport  of  conservative  constituents,  such  as  dye  clouds.  The  unique  as¬ 
pect  of  AdH  is  its  ability  to  ‘adapt’  the  model  domain  mesh  in  certain  loca¬ 
tions  or  at  certain  times  where  more  resolution  may  be  needed.  Sediment 
transport  calculations  are  performed  within  the  SEDLIB  of  AdH.  For  each 
step  of  the  user-specified  time  increment,  AdH  resolves  the  hydrodynam¬ 
ics  (i.e.,  velocities,  water  depths,  etc.)  first,  and  then  passes  that  infor¬ 
mation  to  the  sediment  module.  The  sediment  time-step  is  typically 
smaller  than  the  hydrodynamic  time-step,  i.e.,  multiple  sediment  time- 
steps  are  taken  within  each  hydrodynamic  time-step.  Within  the  sediment 
module,  the  suspended  sediment  concentration  and  bed-load  transport  are 
iteratively  solved  until  they  converge  to  a  change  in  displacement,  which  is 
less  than  the  smallest  grain  diameter.  Once  converged,  the  water  depth  is 
adjusted  based  on  the  converged  displacement  in  order  to  keep  a  constant 
water  surface  elevation,  and  the  water  velocity  is  adjusted  in  order  to  keep 
a  constant  flow  rate.  The  sediment  routine  is  repeated  until  it  reaches  the 
hydrodynamic  time  step,  and  then  the  entire  procedure  is  repeated,  start¬ 
ing  with  the  next  hydrodynamic  time-step. 
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HEC-RAS  is  a  widely  applied  one-  and  two-dimensional  (2D)  unsteady 
flow  model  available  to  the  public  (U.S.  Army  Corps  of  Engineers  2016). 
HEC-RAS  calculates  unsteady  open  channel  hydraulics  using  user-defined 
cross-sections.  For  each  cross-section,  the  one-dimensional  framework 
uses  depth-averaged  and  width-averaged  parameters.  For  each  time-step, 
HEC-RAS  uses  the  sediment  continuity  equation  and  the  sediment 
transport  capacity  to  determine  the  change  in  sediment  volume.  Bed  eleva¬ 
tion  change  and  grain  size  distributions  are  subsequently  calculated  at  all 
nodes  of  each  cross-section. 

1.2  Objective 

The  objective  of  this  literature  review  is  to  prepare  the  conceptual  founda¬ 
tion  for  the  continued  analysis  of  the  flume  experiment  data  and  future 
numerical  modeling  improvements  related  to  bed-load  dispersion.  This  re¬ 
port  provides  a  background  summary  of  dispersion  related  to  solute 
transport,  because  bed-load  dispersion  is  influenced  by  the  same  processes 
as  solute  dispersion.  However,  bed-load  dispersion  is  also  caused  by  addi¬ 
tional  factors  such  as  bed  mixing,  bed  forms,  etc.  which  are  discussed  in 
more  detail  in  Section  3. 

The  report  is  organized  as  follows: 

•  Section  2  explains  the  background  of  diffusion  and  dispersion  con¬ 
cepts. 

•  Section  3  is  focused  on  literature  specifically  related  to  bed-load 
dispersion. 

o  Section  3.1  summarizes  the  available  literature  on  longitudi¬ 
nal  dispersion  using  normal  (Section  3.1.1),  anomalous  (Sec¬ 
tion  3.1.2),  and  fractional  advection-dispersion  (AD) 
relationships. 

o  Section  3.2  provides  a  summary  of  transverse  dispersion  lit¬ 
erature. 
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2  Background 

When  the  applied  skin  friction  bed  shear  stress  exceeds  a  critical  thresh¬ 
old,  non-cohesive  sediment  particles  will  begin  to  move.  Observations  of 
the  transition  show  long  streaks  of  particles  in  motion  very  close  to  the  bed 
with  occasional  bursts  into  the  water  column.  Some  particles  will  be  sus¬ 
pended  into  the  water  column  (suspended  load)  and  others  will  roll,  jump, 
and/or  slide  along  the  bed  (bed-load). 

Sediment  transport  has  been  studied  thoroughly  since  the  seminal  work  of 
Paul  Francois  DuBoys  in  1879.  Early  studies  of  sediment  transport  focused 
on  predicting  the  equilibrium  sediment  transport  rate  and  computing  the 
bed  change  by  solving  a  continuity  equation  (e.g.,  Meyer-Peter  and  Muller, 
1948;  Bagnold,  1956).  These  formulations  computed  mass  change,  but  did 
not  explicitly  account  for  the  stochastic  behavior  of  particles.  Recent  re¬ 
search  suggests  that  the  non-equilibrium  and  stochastic  behavior  of  bed¬ 
load  transport  can  be  predicted  using  AD  models. 

2.1  Diffusion  and  dispersion  concepts 

Diffusion  is  the  process  by  which  a  constituent  ‘spreads’  due  to  molecular 
movements,  collisions,  and  turbulent  fluctuations.  Brown  (1828)  observed 
the  apparently  random  movement  of  pollen  particles  suspended  in  water. 
The  phrase  ‘random  walk’,  first  introduced  by  Karl  Pearson  (1905),  is 
sometimes  used  synonymously  with  Brownian  motion,  where  a  particle’s 
path  consists  of  step  angles  and  distances,  which  are  completely  independ¬ 
ent  from  previous  steps.  Albert  Einstein  was  the  first  to  successfully  ex¬ 
plain  the  cause  of  Brown’s  observations  due  to  the  stochastic  behavior  of 
the  surrounding  fluid  molecules  (Einstein  1905).  Fick  (1855)  identified 
that  in  the  process  of  solute  diffusion  within  a  fluid,  the  flux  is  propor¬ 
tional  to  the  concentration  gradient  of  the  solute,  which  became  known  as 
Fick’s  first  law  of  diffusion  (in  one-dimension): 


(1) 
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where: 

q  =  diffusive  mass  flux  [M  T1  L_1] 

C  =  concentration  [M/L3] 

D  =  diffusion  coefficient  [L2/T] 

‘Dispersion’  also  describes  constituent  ‘spreading’  processes.  These  words, 
diffusion  and  dispersion,  are  often  used  interchangeably  or,  at  least,  with 
overlapping  semantic  range.  However,  dispersion  usually  refers  to  macro¬ 
scale  processes,  which  occur  due  to  the  non-uniform  movements  of  the 
surrounding  fluid.  Holley  (1969)  clarified  the  terminology,  defining  ‘dis¬ 
persion’  as  a  deviation  from  the  temporal  average  convection  velocity,  and 
diffusion  as  the  “difference  between  the  true  convection  in  a  given  direc¬ 
tion  and  the  spatial  average  of  convection  in  that  direction.”  Any  other  lo¬ 
calized  changes  to  the  flow,  such  as  around  structures  or  geometric 
features  of  cross-sections,  can  also  impact  the  dispersion  rate.  Although 
dispersion  can  behave  very  similar  to  diffusion,  dispersive  spreading  is 
generally  one  or  two  orders  of  magnitude  greater  than  molecular  diffusion. 
Additionally,  in  sediment  transport,  where  the  ‘discrete  packets’  of  trans¬ 
portable  mass  (sediment  particles)  are  larger  than  the  molecular  scale, 
typically  only  macro-scale  processes  are  applicable. 

Fick’s  law  of  diffusion  may  accurately  describe  many  scenarios  of  disper¬ 
sion,  but  it  is  generally  an  extrapolation  of  the  original  derivation  based  on 
the  molecule-scale. 

Figure  1  visually  shows  a  simple  numerical  representation  of  how  a  contin¬ 
uous  point  source  could  theoretically  spread  within  a  2-D  rectangular 
channel  having  a  uniform  velocity  from  left  to  right. 
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Figure  1.  Example  of  2-Dimensional  advection-dispersion 
(advection  from  left  to  right). 


2.2  Dispersion  coefficient 

Dispersion  coefficients  (K)  are  difficult  to  measure  or  estimate.  Sir  Geoffrey 
Taylor  was  responsible  for  most  of  the  early  literature  on  these  coefficients, 
including: 

•  Dispersion  by  continuous  movements  (Taylor  1921), 

•  Dispersion  in  laminar  flow  through  a  tube  (Taylor  1953),  and 

•  Dispersion  in  turbulent  flow  through  a  pipe  (Taylor  1954). 

In  general,  Taylor’s  work  explained  how  dispersion  occurs  via  the  velocity 
variations  within  the  fluid  flow.  Elder  (1959)  applied  Taylor’s  analysis  to 
solute  dispersion  within  open-channels  assuming  a  von  Karman  logarith¬ 
mic  velocity  profile  in  order  to  obtain  a  theoretical  dispersion  coefficient 
using  basic  properties  of  the  flow: 

K  =  5.93ut,h  ,  0  , 


where: 

h  =  water  depth  [L] 

u*  =  friction  (or  shear)  velocity  [L/T];  u*  =  y/ghS 
g  =  gravitational  acceleration  [L/T2] 

S  =  channel  slope  [dimensionless] 
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Fischer  et  al.  (1979)  demonstrated  that  dispersion  coefficients  vary  widely 
in  real  streams,  but  are  usually  much  greater  than  Elder’s  result.  For  real 
streams,  Fischer  et  al.  (1979)  developed  the  following  relationship  for  the 
dispersion  coefficient: 


£  =  0.011 


U2W2 

u.,h 


(3) 


where: 


U  =  mean  depth-averaged  velocity  [L/T] 

W  =  channel  width  [L] 

Equation  3  agreed  with  most  study  observations  within  approximately  a  fac¬ 
tor  of  four.  Fisher  et  al.  (1979)  noted  that  there  is  uncertainty  in  determining 
a  specific  dispersion  coefficient  in  real  streams  since  every  irregularity  con¬ 
tributes  to  dispersion,  but  the  majority  of  streams  are  uniform  enough  for 
an  approximate  analysis. 

Many  methods  of  estimating  the  solute  transport  dispersion  coefficient 
have  similar  forms  (McQuivey  and  Keefer  1974;  Liu  1977;  Iwasa  and  Aya 
1991;  Koussis  and  Rodriguez-Mirasol  1998;  Seo  and  Cheong  1998;  Deng 
et  al.  2002;  Kashefipour  and  Falconer  2002;  and  Sahay  and  Dutta  2009), 
with  different  coefficients,  as  shown  in  Table  1.  Disley  et  al.  (2015)  com¬ 
pared  the  equations  listed  in  Table  1  in  more  depth  and  showed  an  im¬ 
proved  behavior  over  a  large  number  of  dye  tracing  field  tests  with  the 
2015  equation  as  listed  in  Table  1. 


Table  1.  Solute  transport  dispersion  coefficient  equations  found  in  the  literature. 


Investigator 

Empirical  Equation  (*) 

Disley  et  al.  (2015) 

K  (  u  \~0A117  /W\0-6776  /u\10132 

uA-3'563Ua)  vJ  U) 

McQuivey  and  Keefer  (1974) 

- (-r  =  0.058  (£)  (—) 
u*h  \SJ  \u */ 

Liu  (1977) 

K  (W\2  fu\os 

u*h~°'18\h)  (uj 

Iwasa  and  Aya  (1991) 

K  =  2  (W) 
u^h  \hJ 
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Investigator 

Empirical  Equation  (*) 

Koussis  and  Rodriguez-Mirasol 
(1998) 

K  (W\2 

u*/i -  0-6  UJ 

Seo  and  Cheong  (1998) 

K  (W\0  62  (  u  \ 1,428 

u,„=5-i95(j  U 

Deng  et  al.  (2002) 

K  (  1  \  /W\ 1,667  /it\0,5 

—  =  0.15  —  —  — 

it*n  V8e/  \  h )  \uj 

(  1  \  /tr\1,38  ( u\ 
e  —  0.145  +  (3^20)  \ji)  [UJ 

Kashefipour  and  Falconer  (2002) 

K  (u\2 

,  =  10.612 

\uj 

Sahayand  Dutta  (2009) 

K  /W\0-96 /u\125 

u^h  \hJ  \uj 

*Note:  See  Appendix  for  symbol  definitions 


The  purpose  of  presenting  the  solute  dispersion  coefficient  equations  in 
Table  l  is  to  compare  the  forms  of  the  equations  that  could  be  useful  for 
understanding  bed-load  dispersion.  Bed-load  dispersion  is  influenced  by 
the  same  factors  as  solute  dispersion,  but  it  is  also  influenced  by  a  number 
of  other  complexities  specific  to  bed-load  sediment  transport,  as  discussed 
in  the  next  section. 
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3  Bed-Load  Dispersion 

The  primary  focus  of  this  report  is  the  dispersion  of  sediment  particles  that 
are  being  transported  along  the  bed.  In  general,  bed-load  dispersion  is 
caused  by  spatial  and  temporal  variations  in  particle  velocities.  The  varia¬ 
tions  in  particle  velocities  can  be  caused  by  many  processes,  including  var¬ 
iations  in  the  advective  and  turbulent  fluid  velocities,  particle  collisions, 
sediment  mixing  in  the  bed,  armoring,  bedforms,  sediment  clusters,  flow 
belts,  channel  morphology,  and  variations  in  particle  characteristics  such 
as  shape,  density,  and  size.  Additionally,  natural  environments  have  many 
other  causes  of  bed-load  dispersion:  vegetation,  structures,  fish,  insects, 
and  macroinvertebrates.  All  of  the  above  processes  cause  bed  sediments  to 
transport  non-uniformly  and  hence  produce  dispersion.  Haschenburger 
(2013)  offers  insight  about  many  causative  factors,  such  as  those  men¬ 
tioned  above,  through  a  17-year  field  study  of  tagged  gravels  in  the  Carna¬ 
tion  Creek. 

Different  modeling  frameworks  may  choose  to  simulate  processes  men¬ 
tioned  in  the  previous  paragraph  independently  from  their  dispersive  bed¬ 
load  transport.  For  example,  simulating  bed  mixing,  sorting,  or  armoring 
separately  can  cause  bed-load  dispersion  within  the  model  on  its  own,  but 
there  is  also  bed-load  dispersion  caused  by  many  processes  that  may  not 
be  modeled,  such  as  listed  in  the  previous  paragraph.  A  representative  dis¬ 
persion  coefficient  that  works  within  one  modeling  framework  may  not  be 
comparable  to  a  dispersion  coefficient  within  a  model  that  handles  the 
processes  differently.  For  example,  the  lvalue  for  bed-load  dispersion  in  a 
model  that  simulates  bed  mixing  as  a  separate  process  is  expected  to  be 
different  from  the  K  value  in  a  model  that  does  not  simulate  bed  mixing. 
Along  these  lines,  the  work  of  Pelosi  et  al.  (2016)  studied  two  different  nu¬ 
merical  modeling  frameworks  of  how  to  handle  vertical  mixing  and 
demonstrated  that  the  second  Exner  Based  Master  Equation  (2D  EBME) 
(Pelosi  et  al.  2014),  which  explicitly  included  vertical  movement,  repre¬ 
sented  field  data  better  than  the  active  layer  formulation  EBME-A.  In  gen¬ 
eral,  a  numerical  dispersion  coefficient  may  vary  significantly  between 
models  and  is  relative  to  the  specific  context  of  each  model’s  considera¬ 
tions. 
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3.1  Longitudinal  dispersion 

3.1.1  Normal  dispersion 

The  literature  diverges  between  those  who  consider  the  bed-load  disper¬ 
sion  as  normal  (Fickian)  behavior,  and  those  who  consider  it  anomalous 
(non-Fickian).  This  section  will  discuss  research  on  normal  dispersion 
while  the  next  section  will  discuss  anomalous  dispersion. 

Under  normal  dispersive  behavior,  a  Gaussian  distribution  typically  repre¬ 
sents  particles  spreading,  and  the  Fickian  analogy  from  diffusion  can  be 
expanded  to  dispersion  of  saltating  grains  existing  in  three  states:  motion, 
stationary,  and  buried.  The  observed  mean  and  variance  of  particle  dis¬ 
placements  from  experimental  measurements  at  specific  time  snap-shots 
can  be  used  to  determine  a  dispersion  coefficient.  The  variance  of  particle 
displacements,  V(x),  can  be  calculated  by: 


V(x)  =  (x  —  x)2  —  x2  —  x2 


(4) 


where: 

=  notation  for  representing  a  mean  value,  averaged  spatially 
x  =  mean  particle  displacement,  averaged  over  the  individual 
particles 

x2  =  mean  of  the  squared  particle  displacements,  averaged  over  the 
individual  particles 

V(x)  is  also  referred  to  as  the  mean  squared  displacement  (MSD).  Einstein 
(1905)  showed  that  the  variance  of  the  particle  displacements  can  be  line¬ 
arly  related  to  time  using  the  dispersion  coefficient  when  taking  a  limit  as 
the  number  of  particle  steps  approaches  infinity.  Following  that  reasoning 
and  using  the  assumption  of  normal  Fickian  dispersion,  the  equation  (Tay¬ 
lor  1921)  for  relating  variance  of  particle  displacements  to  the  dispersion 
coefficient,  K,  is: 


V  =  2  Kt 


(5) 


Several  descriptive  studies  measured  bed-load  tracer  dispersion,  in  the  lab 
(Wong  et  al,  2007;  Gibson,  2009;  Gibson  et  al.,  2011)  and  the  field  (Ken¬ 
nedy  and  Kouba,  1970;  Ferguson  et  al,  2002;  Carling  et  al.,  2006; 
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Haschenburger  2013;  Hassan,  2013),  but  did  not  collect  data  that  could  be 
translated  into  coefficients. 

Hubbell  and  Sayre  (1964)  (Sayre  and  Hubbell  1965)  performed  the  first 
bed-load  dispersion  experiments,  using  radioactive  tracers  to  track  sand 
transport  over  the  course  of  several  days.  As  radioactive  methods  fell  out 
of  favor,  however,  subsequent  studies  focused  on  coarser  gravel  particles 
that  are  easier  to  distinguish  visually. 

Drake  et  al.  (1988)  introduced  125  orange-painted  4-8  mm  particles  onto 
the  bed  of  Duck  Creek  near  Pinedale,  Wyoming,  photographing  their  posi¬ 
tions  every  15  sec  for  240  sec.  Since  the  particles  appeared  to  disperse 
more  quickly  in  the  first  part  of  the  experiment,  dispersion  coefficients 
were  estimated  over  the  final  180  sec  of  the  experiment.  Assuming  normal 
behavior,  Drake  et  al.  computed  dispersion  coefficients  of  approximately 
4.6  cm2/s  in  the  longitudinal  direction  and  0.26  cm2/s  in  the  transverse  di¬ 
rection.  However,  the  caption  of  their  plot  (Drake  et  al.  1988,  and  Figure 
9b)  raises  a  concern  about  the  accuracy  of  their  estimates,  because  they 
state  that  their  dispersion  rates  were  based  on  fitting  straight  lines  “by 
eye”  between  Jv(x)  and  t,  instead  of  between  V(x)  and  t.  It  is  possible  that 
the  estimated  values  of  4.6  and  0.26  cm2/s  were  calculated  correctly,  and 
that  it  is  just  their  figure  caption  which  is  misleading. 

Chang  and  Yen  (2002)  assumed  bed-load  sediment  transport  can  be  mod¬ 
eled  with  the  classic  advection-dispersion  equation  (ADE)  in  one-dimen¬ 
sion: 


dC  dC  v  d2C 

dt  dx  dx  ( 6 ) 


where  u  =  mean  bed-load  material  velocity  [L/T],  K  =  bed-load  dispersion 
coefficient  [L2/T],  and  C  is  the  bed-load  material  concentration  with  units 
usually  expressed  as  mass  per  unit  area  or  volume  per  unit  area.  Chang 
and  Yen  (2002)  estimated  the  dispersion  coefficient  through  dimensional 
analysis,  which  included  the  following  dimensionless  quantities: 


—  =  / 
u,d 


( 


Sne,m, 


(s-\)gd  h 


V 


(7) 
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where: 

Soe  =  initial  equilibrium  bed  slope  [dimensionless] 
m  =  loading  ratio  (calculated  as  the  ratio  of  the  total  sediment 

supply  rate  divided  by  the  sediment  transport  rate  of  the  initial 
equilibrium  state)  [dimensionless] 
v  =  kinematic  viscosity  of  the  fluid  [L2/T] 
s  =  ratio  of  sediment  density  ( ps )  to  water  density  (p) 
[dimensionless] 
d  =  sediment  size  [L] 

Chang  and  Yen  (2002)  assumed  that  the  effect  of  the  sedimentation  pa¬ 
rameter,  v2/[(s  -  1  )gd3],  was  negligible  and,  since  the  range  of  d/h  during 
their  experiments  did  not  change  much,  both  of  those  terms  were  elimi¬ 
nated  from  their  analysis.  They  also  performed  nine  flume  experiments 
with  two  uniform  fine  gravel  particle  sizes  (2.95  mm  and  4.36  mm)  and 
used  a  least  mean  square  error  regression  analysis  to  fit  the  data,  compu¬ 
ting  a  dispersion  coefficient  equation: 

—  =  \06-94  S2oeVsl 

u.,d  (  8  ) 

For  the  range  of  experimental  conditions  tested  within  the  Chang  and  Yen 
(2002)  study,  their  equation  yielded  a  wide  range  of  lvalues,  from 
7.8  cm2/s  to  113  cm2/s. 

A  series  of  papers  by  Furbish  et  al.  (2012  a,  b,  and  c)  and  Roseberry  et  al. 
(2012)  describe  the  theory  and  details  of  bed-load  sediment  flux.  Furbish 
et  al.  (2012c)  use  the  description  of  particle  velocity  autocovariance  (Tay¬ 
lor  1921)  to  calculate  longitudinal  diffusivity  values  in  the  range  of  0.3  to 
0.8  cm2/s.  Although  it  was  intended  as  a  conceptual  demonstration,  as  op¬ 
posed  to  a  process  for  calculating  a  specific  value,  one  section  of  Furbish 
et  al.  (2012b)  back  calculated  a  dispersion  coefficient  of  K  =  3.4  cm2/s 
within  the  Fokker-Planck  formulation.  However,  this  result  depends  on 
simulation  of  specific  assumptions  about  the  travel  distances  for  individual 
particle  trajectories  and  a  certain  impact  on  downstream  particle  activity, 
which  are  uncertain  and  system  specific  parameters.  Furbish  et  al.  (2012b) 
also  explain  how  the  mean  velocity  and  diffusivity  are  fundamentally  re¬ 
lated  to  each  other,  and  how  each  are  connected  to  the  particle  activity 
within  the  bed. 
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Table  2  and  Table  3  summarize  available  relationships  and  data,  respec¬ 
tively,  for  the  bed-load  dispersion  coefficient.  For  the  specific  case  of 
slope-dependent  transport,  Furbish  et  al.  (2009)  developed  a  relationship 
for  bed-load  dispersion  based  on  a  number  of  parameters  as  shown  in  Ta¬ 
ble  2.  Interestingly,  one  of  the  key  fundamental  length  scales  was  the  ac¬ 
tive  soil  thickness,  Has.  Even  more  simple  than  the  Elder  (1959)  equation 
(Eq.  2),  Fan  et  al.  (2014)  used  a  diffusivity  parameter  that  was  propor¬ 
tional  to  the  friction  velocity,  it*,  although  it  has  different  units  (L2T_3)  and 
is  not  exactly  described  as  a  dispersion  coefficient.  Heyman  et  al.  (2015) 
also  found  that  streamwise  dispersion  scaled  linearly  with  shear  velocity 
( K  »  3.3 d50it*).  Table  3  lists  the  available  values  for  the  bed-load  disper¬ 
sion  coefficient  found  in  the  literature. 


Table  2.  Bed-load  dispersion  coefficient  formulas  and  relations  available  in  literature. 


Source 

Formula  or  relation  for  K* 

Notes/limitations 

Furbish 
et  al.  2009 

K  =  ahAS(d/2)Na(  1-^-)  cos 2  9 

\  CmJ 

Diffusion-like  coefficient  in 
slope-dependent  transport 

Chang  and 
Yen  (2002) 

=  10  6-94S02e81m-°-87 
utd 

Based  on  dimensional 
analysis 

Fan  et  al. 
(2014) 

Diffusivity  oc  ut 

note:  this  diffusivity  has 
different  units 

Fleyman  et 
al.  (2015) 

K  ~  3.3 cigoix* 

Based  on  steep,  shallow 
flume  with  uniform  gravel 

*Note:  See  Appendix  for  symbol  definitions 


Table  3.  Bed-load  dispersion  coefficient  values  available  in  literature. 


Source 

Formula  or  value  for  K* 

Notes/limitations 

Furbish  et  al. 

K  =  3.4  cm2/s 

Developed  using  a  Fokker- 

(2012b) 

Planck  formulation 

Furbish  et  al. 

K  =  0.3  to  0.8  cm2/s 

Based  on  Taylor  velocity 

(2012c) 

autocorrelation 

Drake  et  al. 

Kx  =  4.6  cm2/s 

Based  on  field 

(1988) 

Ky  =  0.26  cm2/s 

observations,  using  MSD 

Seizilles  et  al. 

Ky«  0.03 Vs d 

Vs  is  Stoke’s  settling 

(2014) 

velocity 

Fan  et  al.  (2014) 

Diffusivity  oc  ut 

note:  this  diffusivity  has 
different  units 

*Note:  See  Appendix  for  symbol  definitions 
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Most  USACE  developed  transport  models  include  computational  disper¬ 
sion  capabilities;  including  AdH  SEDLIB  (Brown  et  al.  2014),  the  Coastal 
Modeling  System  (CMS)  (Demirbilek  and  Rosati  2011;  Lin  et  al.  2011; 
Sanchez  et  al.  2011a,  and  b),  and  HEC-RAS  (USACE  2016;  Gibson  and 
Boyd  2014).  Even  though  the  hydraulic  routing  elements  in  the  hydrologic 
sediment  yield  model  in  the  HEC-Hydrologic  Modeling  System  (HMS)  in¬ 
cluded  the  (Fischer  et  al.  1979)  analytical  dispersion  equation  in  develop¬ 
mental  versions,  HEC-HMS  did  not  release  that  version  due  to  difficulty 
parameterizing  the  coefficient.  However,  according  to  ERDC  scientists, 
only  AdH  includes  a  bed-load  specific  approach,  which  adds  a  user-speci¬ 
fied  scaling  coefficient  to  Equation  2.  Further  research  is  needed  in  deter¬ 
mining  the  validity  of  the  expression  and  the  appropriate  scaling 
coefficient  for  different  applications. 

3.1.2  Evidence  of  anomalous  dispersion 

Because  of  the  simple  analogy  to  well-established  theory,  bed-load  disper¬ 
sion  and  diffusion  are  usually  discussed  as  representing  normal  (Fickian) 
behavior.  However,  growing  literature  questions  this  analogy,  hypothesiz¬ 
ing  that  bed-load  is  transported  according  to  principles  of  anomalous 
(non-Fickian,  non-normal)  dispersion.  Sun  et  al.  (2015)  flume  data  analy¬ 
sis  showed  that  the  classical  ADE  is  incapable  of  explaining  the  bed-load 
sediment  dispersion. 

Normal  dispersion  affects  the  particle  displacement  variance  proportion¬ 
ally  with  time  (Eq.  5).  However,  anomalous  dispersion  produces  particle 
spreading  rates  either  larger  (super-dispersion)  or  smaller  (sub-disper¬ 
sion)  than  the  normal  model.  Mathematically,  the  general  form  of  relating 
variance,  in  the  x-direction,  to  time,  t,  is: 

V(x)cc  t2y  /Q 


Where: 


y  =  scaling  exponent  for  dispersion 

The  exponent,  y,  is  0.5  for  normal  dispersion,  o  <  y  <  0.5  for  sub-disper¬ 
sion,  or  0.5  <  y  for  super-dispersion. 
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In  addition  to  a  different  relationship  for  the  variance  of  the  particle  dis¬ 
placements,  anomalous  dispersion  can  also  affect  the  shape  of  the  proba¬ 
bility  density  function.  Normal  dispersion  typically  approximates  a 
Gaussian  distribution,  as  shown  in  Figure  2.  With  this  assumption,  the 
peak  drops  and  probabilities  widen  as  time  progresses,  but  the  particle 
distribution  remains  symmetrical  around  the  mean.  However,  it  is  possi¬ 
ble  for  other  types  of  distributions  that  are  not  symmetrical  to  describe  the 
dispersion  process.  Indeed,  Grigg  (1969)  and  Yang  and  Sayre  (1971)  found 
that  the  stochastic  behavior  of  sediment  particles  behaved  according  to  a 
skewed  gamma  distribution,  more  like  Figure  3,  as  opposed  to  a  shape 
similar  to  those  in  Figure  2. 

Figure  2.  Example  of  dispersion  that  follows  a  Gaussian  distribution,  while  also 

undergoing  advection  from  left  to  right  (progressing  in  time  from  ti  ->  h  ->  fc). 


Even  during  solute  transport,  anomalous  dispersion  occurs  initially.  Ac¬ 
cording  to  Fischer  et  al.  (1979),  normal  dispersion  does  not  apply  at  nor¬ 
malized  distances  ( x'  —  xet/(uW2) )  less  than  0.4,  where  et  is  a  transverse 
mixing  coefficient.  In  the  initial  spreading  phase,  the  distribution  can 
skew,  pushing  the  peak  left  of  the  mean  and  elongating  the  leading  edge 
(i.e.,  a  heavy  tail)  (Hill  et  al.  2010),  as  sketched  in  Figure  3.  Liebault  et  al. 
(2012),  claims  that  frontrunners,  or  particles  at  the  leading  edge,  can 
travel  very  quickly,  resulting  in  the  heavy-leading-tail  distribution  shape. 
Bed-load  transport  literature  refers  to  several  types  of  heavy-tailed  distri¬ 
butions  such  as  lognormal,  inverse  gamma,  exponential,  power-law, 
Weibull,  Pareto,  and  Poisson  (Einstein  1972;  Hassan  et  al.  1991;  Hill  et  al. 
2010;  Haschenburger  2013).  Hubbell  and  Sayre’s  (1964)  data  also  showed 
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a  heavy-leading-tail  distribution.  Identifying  which  distribution(s)  best  de¬ 
scribe  bed-load  dispersion  is  an  active  area  of  research  and  an  important 
consideration  when  applying  dispersion.  Fischer  et  al.  (1979)  acknowl¬ 
edged  that  the  longitudinal  concentration  distribution  of  solutes  only  be¬ 
comes  normal  asymptotically,  and  for  normalized  distances  greater  than 
x’  =  1.0  it  is  a  sufficient  approximation.  However,  the  work  of  Phillips  et  al. 
(2013)  suggested  that  the  asymptotic  limit  of  bed-load  dispersion  is  super- 
dispersive.  Fischer  et  al.  (1979)  acknowledged  two  examples  of  dispersion 
studies  that  reportedly  never  showed  normal,  Fickian  behavior  (Nordin 
and  Sabol  1974,  and  Day  1975). 

Figure  3.  Example  of  dispersion  that  produces  a  heavy-tail  distribution  (progressing  in 

time  from  ti  ->  ->  fe). 


To  investigate  normal  vs.  anomalous  dispersion,  Hassan  et  al.  (2013)  stud¬ 
ied  64  gravel-bed  field  tracer  experiments  to  determine  whether  the  trac¬ 
ers  exhibited  heavy-tail  or  thin-tail  distributions.  Thin-tail  distributions 
indicate  normal  dispersive  behavior,  while  heavy-tail  distributions  indi¬ 
cate  anomalous  dispersion.  Of  the  64  surveys,  51  exhibited  thin-tail  distri¬ 
butions  and  8  more  could  have  been  considered  thin-tail  based  on  the 
definition  of  the  ‘tail’.  Liebault  et  al.  (2012)  performed  experimental  work 
on  the  Bouinenc  Torrent,  a  tributary  to  the  Bleone  River  in  southeast 
France.  Their  work  classified  two  out  of  three  Bouinenc  Torrent  experi¬ 
mental  cases  as  heavy-tail  distributions,  but  upon  reanalysis,  Hassan  et  al. 
(2013)  classified  all  three  Bouinenc  Torrent  surveys  as  thin-tail  distribu¬ 
tions.  Therefore,  the  work  of  Hassan  et  al.  (2013)  indicates  normal  disper¬ 
sive  behavior  in  most  cases,  but  much  of  the  literature  does  not  agree. 

Recent  non-Fickian  literature  (e.g.,  Schumer  et  al.  2009;  Bradley  et  al. 
2010;  Martin  et  al.  2012)  points  to  the  work  of  Nikora  et  al.  (2001  and 
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2002)  as  a  key  concept-framing  work  for  anomalous  bed-load  dispersion. 
Nikora  et  al.  (2001)  identified  three  different  scales  of  bed-load  transport; 
namely  the  local,  intermediate,  and  global  scales.  The  local  scale  is  an  indi¬ 
vidual  particle’s  trajectory  within  one  saltating  movement.  The  particle 
trajectory  between  two  successive  periods  of  rest  is  the  intermediate  scale. 
The  third  scale,  combining  multiple  intermediate  scales,  is  the  global  scale. 
Nikora  et  al.  (2002)  suggest  that  each  range  can  exhibit  a  different  type  of 
dispersive  behavior,  i.e.,  each  can  have  a  different  exponent,  y,  in  Equation 
9.  Bialik  et  al.  (2015)  estimates  the  transition  points  between  scales. 

3. 1.2.1  Local  scale 

The  dispersive  spreading  within  the  local  range  is  caused  by  the  variations 
of  particle  velocities,  accelerations,  hop  distances,  and  travel  times  within 
a  single  jump  (Fathel  et  al.  2015).  There  are  a  number  of  different  types  of 
distributions  used  to  characterize  the  velocities  and/or  travel  times  at  the 
local  scale,  such  as  an  exponential  distribution  (Roseberry  et  al.  2012,  Fur¬ 
bish  et  al.  2012b,  and  Lajeunesse  et  al.  2010)  or  a  truncated  Gaussian  dis¬ 
tribution  (Martin  et  al.  2012,  Ancey  and  Heyman  2014).  The  local  scale  is 
generally  considered  super-dispersive  (Khzeri  and  Chanson  2015,  Bialik 
et  al.  2015).  Martin  et  al.  (2012)  studied  the  physical  mechanisms  at  the  lo¬ 
cal  scale,  observing  super-dispersion  with  an  approximate  exponent,  y,  of 
0.8  at  times  from  0.2  to  1  sec  and  then  y  =  0.7  beyond  1  sec  in  their  experi¬ 
ments.  These  exponent  values  were  in  good  agreement  with  the  super-dis¬ 
persive  observations  of  Nikora  et  al.  (2002),  where  y  was  0.77  to  0.87  for 
the  Balmoral  Irrigation  Canal  in  North  Canterbury,  New  Zealand.  Martin 
et  al.  (2012)  hypothesized  that  the  super-dispersive  behavior  at  the  local 
scale  is  caused  by  particle  inertia  and  fluid  drag,  along  with  possible  influ¬ 
ences  from  hydrodynamic  and  boundary  effects. 

Huang  et  al.  (2011)  studied  the  displacement  of  particles  at  the  molecule- 
scale  and  observed  the  full  transition  from  pure  ballistic  motion  (y  =  1)  to 
diffusive  Brownian  motion  as  the  movement  became  more  defined  by  the 
fluid  behavior.  However,  a  growing  amount  of  literature  points  to  super- 
ballistic  (y  >  1)  motion  for  bed  sediment  dispersion  (Campagnol  et  al. 

2015,  Olinde  and  Johnson  2015,  Phillips  et  al.  2013,  Martin  et  al.  2012  and 
2014).  Campagnol  et  al.  (2015)  showed  that  if  the  particles  are  exhibiting 
short  saltating  hops  (i.e.,  weak  transport  conditions),  the  super-ballistic 
behavior  at  the  local  scale  is  from  the  strong  acceleration  experienced  by 
the  particles. 


ERDC/CHL  SR-16-5 


18 


3. 1.2.2  Intermediate  scale 

The  intermediate  scale  is  generally  characterized  by  the  particle  wait  times 
between  saltating  movements.  In  the  words  of  H.  Einstein  (1942):  “bed¬ 
load  movement  is  to  be  considered  as  the  motion  of  bed  particles  in  quick 
steps  with  comparatively  long  intermediate  periods  of  rest”  (Nikora  et  al. 
2001).  Particle  wait  times  were  significant  enough  in  the  work  of  Bennett 
and  Nordin  (1977)  that  the  researchers  considered  dispersion  of  sus¬ 
pended  sediment  negligible  with  respect  to  dispersion  caused  by  storage  of 
particles  in  the  streambed. 

Literature  offers  conflicting  explanations  for  dispersive  behavior  in  the  in¬ 
termediate  scale.  Nikora  et  al.  (2002)  stated  that  the  intermediate  scale 
could  exhibit  normal,  super-,  or  sub-dispersive  behavior.  The  researchers’ 
analysis  of  Balmoral  Canal  (North  Canterbury,  New  Zealand)  showed  su¬ 
per-dispersion  in  the  intermediate  scale,  but  Zhang  et  al.  (2012)  later  re¬ 
classified  this  work  as  the  local  scale.  Martin  et  al.  (2012)  explained 
challenges  with  identifying  the  transition  between  local  and  intermediate 
scales,  but  stated  that  the  heavy-tailed  wait  times  cause  the  behavior  to  be¬ 
come  sub-dispersive  beyond  the  local  scale.  Zhang  et  al.  (2012)  also  found 
sub-dispersive  behavior  in  the  intermediate  range.  Zhang  et  al.  (2012) 
point  out  that  experimental  observations  can  sometimes  be  challenging 
since  newly  introduced  tracer  particles  tend  to  be  more  mobile  than  aver¬ 
age,  causing  super-dispersion  to  be  overestimated,  while  on  the  other 
hand,  the  transport  of  the  most  mobile  particles  can  be  limited  by  the  win¬ 
dow  of  observation,  which  could  cause  super-dispersion  to  be  underesti¬ 
mated.  Perhaps  this  helps  explain  some  of  the  conflicting  evidence  found 
in  the  literature. 

It  is  important  to  note  that  the  vertical  exchange  within  the  bed  has  an  im¬ 
portant  influence  on  the  particle  storage  and  release,  or  the  particle  wait 
time  between  jumps  (Parker  et  al.  2000;  Pelosi  et  al.  2016).  Dispersion  at 
intermediate  and  global  scales  depends  on  this  vertical  storage  and  ex¬ 
change.  Therefore,  a  sediment  model  could  completely  ignore  the  disper¬ 
sion  term  of  the  ADE,  but  if  it  simulates  vertical  exchange  behavior,  the 
model  results  will  still  exhibit  a  dispersive  effect  (Ferguson  and  Hoey 
2002;  Ganti  et  al.  2010). 
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3. 1.2.3  Global  scale 

The  global  scale  is  the  largest  of  the  three  spatial  scales,  and  it  is  defined  as 
the  combination  of  many  intermediate  scales.  The  sub-,  normal,  or  super¬ 
dispersion  of  the  global  range  is  even  more  controversial  than  agreement 
at  the  intermediate  scale.  Meerschaert  et  al.  (2008)  proposed  a  tempered 
model  to  capture  the  slow  transition  from  sub-dispersion  to  a  normal  dis¬ 
persion  limit.  Voepel  et  al.  (2013)  studied  waiting  times  and  the  transition 
between  the  intermediate  and  global  scales,  reporting  transitions  on  the 
order  of  seconds-to-minutes  in  flume  studies,  while  it  may  take  days  to 
years  in  natural  channels.  Drake  et  al.  (1988),  assuming  normal  disper¬ 
sion,  was  the  first  to  analyze  data  from  Duck  Creek  near  Pinedale,  Wyo¬ 
ming.  Then,  Nikora  et  al.  (2002)  reanalyzed  the  data  to  find  sub-dispersive 
behavior  and  classified  it  in  the  global  range.  Zhang  et  al.  (2012)  advanced 
Meerschaert  et  al.’s  (2008)  tempered  stable  model  and  reanalyzed  the 
Duck  Creek  data  further.  Zhang  et  al.  (2012)  agreed  with  Nikora  et  al. 
(2002)  that  it  was  sub-dispersive,  but  re-classified  the  range  as  intermedi¬ 
ate  instead  of  global.  Bialik  et  al.  (2012)  confirmed  the  existence  of  at  least 
three  ranges  of  scales  and  identified  dispersive  exponents  that  decreased 
from  one  range  to  the  next.  However,  Zhang  et  al.  (2012)  might  classify  all 
of  Bialik  et  al.’s  (2012)  observations  within  the  local  and  intermediate 
ranges. 

Sayre  and  Hubbell  (1965)  performed  field  experiments  using  40  pounds  of 
radioactive  tracers  in  the  North  Loup  River,  Nebraska.  Their  observations 
revealed  an  episodic  behavior  in  which  tracer  particles  were  buried  under 
migrating  dunes  and  then  released  by  scour.  Assuming  exponential  distri¬ 
butions  for  step  distance  and  rest  periods,  Sayre  and  Hubbell  (1965),  in 
line  with  Einstein  (1972),  derived  concentration  profiles.  However,  the  ob¬ 
served  plumes  had  heavier  tails  than  the  Sayre  and  Hubbell  model  indi¬ 
cated,  and  Bradley  et  al.  (2010)  and  Zhang  et  al.  (2012)  proposed  new 
models  in  an  attempt  to  simulate  the  Sayre  and  Hubbell  (1965)  observa¬ 
tions. 

Olinde  and  Johnson  (2015)  used  radio  frequency  identification  and  accel¬ 
erometer  embedded  tracers  to  show  that  for  mountain  streams,  the  vari¬ 
ance  can  actually  begin  to  decrease  after  some  amount  of  time  due  to 
spatial  changes  in  grain  sizes  and  reach  morphology.  Olinde  and  Johnson 
(2015)  found  that  the  displacements  and  step  lengths  exhibited  thin-tailed 
distributions,  while  rest  times  in  mountainous,  snow-impacted  streams 
were  described  by  heavy-tailed  power  law  distributions  resulting  in  an 
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overall  super-dispersive  behavior  of  the  tracers.  Pelosi  et  al.  (2016)  showed 
that  the  2D  EBME  model,  which  resembled  field  data  well,  also  exhibited 
super-dispersive  behavior.  Therefore,  there  is  literature  evidence  for  sub-, 
normal,  and  super-dispersion  at  the  global  scale. 

3.1.3  Fractional  Advection-Dispersion  Equation  (ADE) 

Recent  bed-load  dispersion  modeling  literature  investigates  the  difference 
between  anomalous  dispersion  and  the  classic  ADE.  Specifically,  the  dis¬ 
persion  term  of  the  ADE  for  anomalous  dispersive  behavior  has  a  frac¬ 
tional  derivative,  yielding  the  fractional  advection  dispersion  equation 
(FADE): 


dC  dC  d(l,r)C 

- b  U - —  K - 777 

dt  dx  dx0  y  (  10 ) 

where  y  is  the  scaling  exponent,  the  dispersion  coefficient  K  has  dimen¬ 
sions  of  [LP/r)/T],  and  u  is  the  mean  bed-load  material  velocity.  As  before, 
y  =  0.5  yields  the  normal,  or  Fickian,  dispersive  behavior  and  the  classic 
ADE.  The  work  of  Sabzikar  et  al.  (2015)  applied  tempered  fractional  calcu¬ 
lus  to  describe  the  complexities  of  the  mathematical  behavior.  Benson 
et  al.  (2013)  connected  the  FADE  to  random  walks  governed  by  heavy¬ 
tailed  probability  distributions  and  discussed  several  Eulerian  and  Lagran- 
gian  numerical  methods  for  treating  the  fractional  partial  derivative. 

Ganti  et  al.  (2010)  demonstrated  how  a  formulation  of  the  Exner  equation 
for  sediment  conservation,  which  incorporates  a  probability  distribution  of 
step  lengths,  having  heavy-tails  in  real  rivers,  leads  to  an  anomalous  ad- 
vection-dispersion  formulation.  Fan  et  al.  (2014)  developed  a  Langevin 
formulation  of  bed-load  modeling  using  a  balance  of  the  forces  on  parti¬ 
cles.  Pelosi  and  Parker  (2014)  offered  an  informative  comparison  between 
numerical  frameworks  based  on  flux  and  those  based  on  entrainment. 
Pelosi  et  al.  (2014  and  2016)  used  the  Parker  et  al.  (2000)  framework  of 
the  Exner  equation  to  simulate  complexities  within  the  bed,  departing 
from  the  ‘active  layer’  concept  and  producing  vertical  dispersion. 

Sun  et  al.  (2015)  also  showed  that  classical  ADE  was  incapable  of  explain¬ 
ing  the  behavior  of  laboratory  flume  dispersion,  and  the  researchers’  anal¬ 
ysis  used  fractional  derivatives  for  the  time  derivative  (left-most  term  of 
Eq.  10)  as  well  as  the  spatial  derivative  (right-most  term  of  Eq.  10).  Once 
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implemented,  the  Sun  et  al.  (2014)  FADE  model  was  relatively  insensitive 
to  the  dispersion  coefficient  chosen;  the  tested  range  covered  two  orders  of 
magnitude.  Zhang  et  al.  (2014)  developed  the  use  of  a  simpler  (than 
FADE)  approach  called  the  subordinated  advection  equation  (SAE)  to  ex¬ 
plain  various  historical  data  sets  for  sediment  dispersion  across  a  wide 
range  of  scales.  Numerous  other  recent  numerical  models  targeted  sto¬ 
chastic,  anomalous  dispersion  of  bed-load  transport  including:  Schumer 
et  al.  (2009),  Bradley  et  al.  (2010),  Lajeunesse  et  al.  (2010),  Zhang  et  al. 
(2012),  Lajeunesse  et  al.  (2013),  Martin  et  al.  (2014),  Ancey  and  Heyman 
(2014),  and  others. 

3.2  Transverse  dispersion 

If  the  dispersion  behavior  is  the  same  in  the  longitudinal  (i.e.,  streamwise) 
and  transverse  (i.e.,  cross-stream)  directions,  it  is  considered  isotropic; 
otherwise  it  is  anisotropic.  The  majority  of  the  recent  literature  considers 
it  anisotropic  (Furbish  et  al.  2012c;  Nikora  et  al.  2001;  Roseberry  et  al. 
2012;  Seizilles  et  al.  2014;  etc.).  Samson  et  al.  (1998)  studied  the  motion  of 
spheres  down  a  sloped  bumpy  plane  and  found  that  the  dispersion  could 
be  described  by  anisotropic  normal  dispersion.  Samson  et  al.  (1998)  found 
that  the  longitudinal  dispersion  coefficient  had  a  constant  value  independ¬ 
ent  of  the  sphere  size  and  slope,  while  the  transverse  dispersion  coefficient 
was  found  to  decrease  with  the  plane  slope  and  scale  with  the  size  of  the 
spheres.  Interestingly,  Samson  et  al.  (1998)  found  that  the  longitudinal 
dispersion  coefficient  observed  was  smaller  than  the  transverse  coefficient, 
a  finding  contradicted  by  the  results  of  Roseberry  et  al.  (2012)  where  the 
longitudinal  dispersion  coefficient  was  found  to  be  larger. 

Seizilles  et  al.  (2014)  also  found  that  cross-stream  dispersion  was  entirely 
consistent  with  normal  Fickian  behavior,  in  contrast  to  the  observations  of 
Nikora  et  al.  (2002).  Bialik  et  al.  (2012)  found  only  two  scales  of  disper¬ 
sion  in  the  transverse  direction,  compared  to  at  least  three  in  the  longitu¬ 
dinal  direction.  Similar  to  Samson  et  al.  (1998),  Bialik  et  al.  (2012)  stated 
that  the  spread  of  particles  in  the  longitudinal  direction  mostly  depended 
on  turbulence,  instead  of  relative  particle  size,  while  spreading  in  the 
transverse  direction  did  depend  on  relative  particle  size.  Expanding  upon 
that,  Barati  et  al.  (2015)  used  numerical  modeling  to  show  that  the  main 
mechanism  for  longitudinal  dispersion  is  the  flow  condition,  while  the 
main  mechanisms  for  transverse  dispersion  are  collisions,  both  inter-par¬ 
ticle  and  particle-bed. 
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4  Conclusion 

Many  investigators  have  successfully  applied  Fickian  models  to  their  data, 
but  growing  evidence  suggests  that  non-Fickian,  anomalous,  anisotropic 
processes  offer  value  in  representing  bed-load  dispersion.  Some  research¬ 
ers  agree  that  dispersion  can  be  characterized  by  three  different  scales;  lo¬ 
cal,  intermediate,  and  global.  A  majority  agree  that  the  local  scale  is 
characterized  by  super-dispersive  behavior.  There  are  numerous  findings, 
even  in  just  the  last  couple  of  years,  which  present  significant  disagree¬ 
ments  about  the  global  scale,  whether  it  exhibits  super-,  sub-,  or  normal 
dispersion.  As  witnessed  by  the  number  of  recent  publications,  sediment 
bed-load  dispersion  is  an  active  area  of  research,  and  the  published  values 
for  the  bed-load  dispersion  coefficient  vary  greatly.  Some  formulas  and 
recommended  values  for  the  dispersion  coefficient  are  shown  in  Table  2 
and  Table  3.  For  numerical  modeling,  numerical  dispersion  coefficients 
are  relative  to  the  context  of  all  considerations  within  a  modeling  frame¬ 
work  and  may  not  be  directly  comparable  across  studies.  Future  work  in¬ 
cludes  the  evaluation  and  improvement  of  USACE  sediment  transport 
models  using  the  laboratory  experiments  of  Gibson  et  al.  (2016)  to  test  the 
dispersion  frameworks  discussed  herein. 
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Appendix:  Definitions  of  Symbols 


Symbol  Description  [units;  M=mass,  L  =  length ,  T=time] 


A 

C 

c 


d 

dso 

D 

9 

h 

Has 

K 

K, 

Ky 

m 

MSD 

Na 

9 

S 

s 

See 

At 

£i>  £2 

U 

u 

u* 

V(x) 

Vs 

w 

x 

x’ 

(x) 

a 

Y 

Yx,  Yy 

£t 

e 

K 

V 

P 

Ps 


Area  of  control  volume  normal  to  the  flux  [L2] 

Concentration  (sediment  or  solute)  [typically  M  /  L3] 

Volumetric  particle  concentration  [dimensionless] 

Natural  consolidated  (maximum)  volumetric  particle  concentration 
(sediment  or  solute)  [dimensionless] 

Sediment  size  [L] 

Median  grain  diameter  [L] 

Diffusion  coefficient  [typically  L2  /  T] 
gravitational  acceleration  [L  /  T2] 

Fluid  depth  [L] 

Active  soil  layer  thickness  [L] 

Dispersion  coefficient  [typically  L2  /  T] 

Dispersion  coefficient  in  the  x-direction  [typically  L2  /  T] 
Dispersion  coefficient  in  the  y-direction  [typically  L2  /  T] 

Loading  ratio  of  total  sediment  supply  rate  to  the  sediment 
transport  rate  of  the  initial  equilibrium  state  [dimensionless] 

Mean  Squared  Displacement  [L2] 

Particle  activation  rate  [T1] 

Solute  mass  flux  [M/(L2T)] 

Slope  of  the  channel  [dimensionless] 

Ratio  of  sediment  density  to  water  density  [dimensionless] 

Initial  equilibrium  bed  slope  [dimensionless] 

Time  step  [T] 

Mean  times  of  passage  of  tracer  at  locations  xi  and  x2  [T] 
Stream-wise  velocity  [L  /  T] 

Mean  stream-wise  velocity  [L  /  T] 

Friction  velocity  [L  /  T] 

Variance  of  particle  displacements  in  the  x-direction  [L2] 

Stoke’s  settling  velocity  [L  /  T] 

Width  of  the  channel  [L] 

Particle  displacement  [L] 

Normalized  distance  [dimensionless] 

Mean  particle  displacement  [L] 

Coefficient  [dimensionless] 

Scaling  exponent  for  diffusive  behavior  [dimensionless] 

Diffusive  exponents  in  the  x-  and  y-directions  [dimensionless] 
Transverse  mixing  coefficient  [L2/  T] 

Angle  of  land-surface  slope  [dimensionless] 

Diffusivity  [L2  /  T] 

Kinematic  viscosity  of  the  fluid  [L2  /  T] 

Fluid  density  [M  /  L3] 

Sediment  density  [M  /  L3] 
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